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ABSTRACT 

In a microlensing event, a large magnification occurs at caustic crossing and provides an 
opportunity to obtain a stronger signal associated with the object. In this paper we study the 
possibility of magnetic field detection in a microlensing event through the Zeeman effect. 
We follow the prescription introduced by (Robinson 1980) which analyses the spectrum of a 
star in a Fourier space to deconvolve other broadening mechanism from the actual Zeeman 
effect. First we study magnification contrast between source and spot in terms of spot size and 
then consider two distinct strategies using modern spectrographs and perform a Monte Carlo 
simulation to find the detection efficiency in magnification-magnetic field plan. The spectral 
resolution and signal to noise ratio in each strategies, specify suitable places in this plane 
to detect a magnetic field. Apart from complexity of magnetic field on a star, we consider 
a simple model and propose a fantastic method to detect magnetic field using spectrum of 
source stars at caustic crossing. 

Key words: Physical Data and Processes: gravitational lensing: micro - magnetic fields, 
Stars: activity, Astronomical techniques: spectroscopic 


1 INTRODUCTION 

When a compact object (lens) sufficiently aligned with the light 
source, two images of the source appear due to light bending in 
the gravitational field. For the galactic events the angular distance 
of images is of the order of mili-arc second, thus can’t be resolved 
using ground based instruments and called micolensing. This mul¬ 
tiplication of images manifests itself as a magnification of the 
source flux and provides a powerful tool to study faint objects. 
For a point like lens the magnification of a point like source is 
an analytic function but magnification of a uniformly bright cir¬ 
cular source is given by elliptical integrals(Witt & Mao 1994). In 
contrast to a point lens, magnification for a binary system can’t 
be expressed by analytic function and two main approaches have 
been developed to find the magnification, ray-shooting algorithm 
(Kayser et al. 1986a; Schneider & Weiss 1987) and adoptive coun¬ 
tering (Dominik 2007). 

Gravitational microlensing provides a unique tool to study 
not only the properties of the lens system, (e.g extra-solar planets 
(Cottle et al. 1991; Gould & Loeb 1992; Gould 2009)), but also the 
properties of the source star (e.g stellar atmosphere (Albrow et al. 
1999; Afonso 2001; Gould 2001; Abe et al. 2003; Fields et al. 
2003)). Recently author of (Rahvar 2015) has collected the most 
recent application of microlensing. For a binary or multiple point 
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like lens systems there are some closed curved in the source plane 
which the magnification of a point source diverges and called caus¬ 
tic line. The corresponding images position is called critical curved 
and formed closed curve like caustic lines. (For more information 
see (Erdl & Schneider 1993) ). Near a caustic line an extended 
source undergoes a large magnification and due to high magnifi¬ 
cation gradient, the source surface is differentially magnified. Rel¬ 
atively slow lens-source proper motion, provides an opportunity 
to sample the light curve with high cadence and study detail of a 
source surface. For example the technique has been used to obtain 
limb-darkening profile of Bulge giant and sub-giant (Fields et al. 
2003) as well as main sequences (Abe et al. 2003). 

In this paper we study the possibility of measuring the mag¬ 
netic field which may be existed at the surface of a source star, 
when the source is crossing the caustic line. Caustic crossing has 
two major impacts on the embedded magnetic field, on one hand 
the spot area will increase and on the other hand total flux is mag¬ 
nified. These two effects increase the signal to noise ratio (SNR) 
and provide a special time to detect the magnetic field. In the next 
paragraph we will shortly discuss the detection of stellar magnetic 
field. 

Evidences for the presence of the Solar magnetic field was 
established by Hale in 1908 (Hale 1908). He measured the mag¬ 
netic filed of sunspots to be ~ 2600 — 2900 Gauss using the line 
splitting and polarization. Detection of the magnetic field over the 
surface of other stars is much more difficult. There are several phe- 
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nomenons in the range of magnetic activity that might mimic mag¬ 
netic field. These effects include dark spots (Berdyugina 2005), 
flares (Favata et al. 2000), enhanced chromospheric emission in 
some bands (Pizzolato et al. 2003; Hall & Lockwood 2004), and 
hence, distinguishing magnetic field from these effects is a difficult 
task. In spite of such difficulty, average magnetic field of Ap stars 
was measured using Zeeman broadening in 1971 (Preston 1971). 
(Robinson 1980) developed an interesting technique to measure 
the amplitude of the magnetic field using Fourier analysis of the 
line profile produced due to Zeeman effect. His method involved 
crude assumptions regarding the radiative transfer which later im¬ 
proved by other authors. (Saar 1988; Basri & Marcy 1988; Saar 
1990; Ruedi et al. 1997). The detection or measuring the magnetic 
field has been revolutionised over the past decades due to the de¬ 
velopments of the new generation of spectrographs (VLT/UVES, 
ESO/HARPS, Keck/HIRES, etc.) and spectropolarimeters Narval 
at TBL and HARPSpol at ESO. More advanced discussions regard¬ 
ing this issue can be found in (Mathys 2012) and (Landstreet 2013). 
To best of our knowledge, there is no attempt to find magnetic field 
in a microlensing event. In this work apart from all complexity, we 
use a simple model for magnetic field on the source star and study 
the possibility of its detection when the source is crossing a caustic 
line. 

The structure of this paper is as follows. In Sec.(2) we present 
the key equation of microlensing and amplification at caustic line. 
Then we investigate magnification contrast between magnetic and 
non magnetic area and perform a Monte Carlo simulation to show 
how much the area of a spot will increase at caustic crossing. In 
Sec.(3) we show how a spectral line be modified in the presence of 
a magnetic field and how one can detect the signal of the magnetic 
field in the spectrum of a star. In Sec.(4) the possibility of mag¬ 
netic field detection is studied using Monte Carlo simulation and in 
Sec.(5) we conclude and discuss the results. 


2 GRAVITATIONAL MICROLENSING AND CAUSTIC 
LINE 


Gravitational field of a lens star deflects light of a source star. Such 
a deflection, depending on the mass distribution of lens, may pro¬ 
duces two or several images. Since the gravitational field doesn’t 
change the surface brightness of the source, magnification is sim¬ 
ply the ratio of images to source area. 

Using thin sheet approximation for the lens, the deflection an¬ 
gle is given by (Schneider 1985): 



where x is a 2 dimensional vector normalized to 
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which is the Einstein radius in the lens plane. Where D s ( Da ) are 
the source (lens) distance respectively and p — In Eq. (1) the 
definition of parameters is as following, 
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where £(as) is the lens surface mass density. The lens equation 
which describe the location of images is given by: 


y = x - a(x) , 


(4) 


where y shows a point source in the source plane and normalized 
to r e . For a binary system, the lens equation is: 
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where x± and X 2 are the location of lenses. If the origin of coordi¬ 
nate system set to the center of mass we have 
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Where d and q are the distance and mass ratio of lenses and with¬ 
out loss of generality we put the lenses along horizontal axes. In 
addition to d and q parameters, we need three parameters includ¬ 
ing uo (minimum impact parameter), to (time to reach uo) and ts 
(The Einstein crossing time) to completely describe a binary event. 
The magnification of each image is simply the ratio of image to the 
source area and is given by 

Pi = det(^)~ 1 (x i ) , (7) 

and the total magnification is a sum over individual magnification, 

p(y) = ) _ 1 (*i) • (8) 


Critical points is defined via 

d e t(^)(x c ) = 0 , (9) 

and form closed curves. On the other hand all points in the source 
plane y c = y(x c ), also form closed curves which so called caustic 
line. When a point source crosses a caustic, according to Eq.(7), 
its magnification diverges, considering finite size effect the mag¬ 
nification doesn’t diverge but become very large. The magnifica¬ 
tion of an extended source can be obtained via inverse rate shoot¬ 
ing algorithm(Kayser et al. 1986b). However Dominik introduced 
an efficient algorithm for calculating the magnification(Dominik 
2007). we use (Dominik 2007) algorithm to find the magnification 
throughout of this paper. 

Now we consider a simple model to study the magnetic field 
detection. We consider a source star with a spot which cover certain 
amount of projected source area. We parametrize each spot with 
f — which show the initial coverage of spot before crossing 
the caustic, where A s (A u ) is the area of spot (source). 

When such sources cross the caustic line, all parts of the 
source magnifies but a magnification contrast occur between re¬ 
gions close and far from caustic. Such a magnification contrast in¬ 
creases the signal of the magnetic field when the spot is over caus¬ 
tic so makes it more detectable. We consider a source with a sin¬ 
gle spot and define the contrast as —. Where li s is the magnifica- 
tion of the spotted area and p u is the magnification of a uniform 
source. In Fig. 1 the magnification contrast is shown for 4 spots 
with different sizes. As we expected, decreasing the spot size in¬ 
creases the magnification contrast but the initial signal is less than 
a large spot. The amplitude of magnification contrast depends on 
the location of the spot at caustic crossing, we perform a Monte- 
Carlo simulation to find probability of occurring a specific magni¬ 
fication contrast for spots with different sizes. In our analyzes, we 
simulate several events with a spot which located over the source 
randomly and the magnification contrast is average of this quantity 
in each event. In this analyze we select a typical binary lenses with 
(d = 1, q — 0.001, tE = 25 days) and as long as we are near 
central caustic, our results are not sensitive to change the param¬ 
eters of lenses. The results of such an exercise for four different 
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Figure 1. Magnification contrast for spots in 4 different sizes in a typical 
binary lenses event. The red solid, dashed green, tiny dashed blue and dotted 
pink curves represent magnification contrasts for spot with / = 0.001, 
/ = 0.005, / = 0.01 and f = 0.05 respectively. 



Figure 2. Frequency of contrast for spots with different sizes. The mean 
contrast is calculated as a ratio of spotted to the uniform source magnifica¬ 
tion around the caustic line (see text for more detail). The contrast for small 
spot might be larger than a big spot. 


/ is presented in Fig(2). The contrast can be as large as 10 for a 
spot with initial size 0.001 while for a large spot high contrasts are 
less likely. In addition to this effect, when a star crosses a caustic 
line, its brightness might magnify several hundreds times and so 
the signal will increase significantly. Following, after introducing 
two strategies, we consider this effect to detect the magnetic field. 


3 SPECTRAL LINE DISTORTION IN THE PRESENCE 
OF MAGNETIC FIELD 

Strong magnetic field could be produced at the surface of stars due 
to rotation of interior plasma. These magnetic fields split a sensitive 
line due to Zeeman effect. A spectral emission or absorption line 
is excited if electron makes a transition from one state to another. 
Electron in each energy state has an orbital angular momentum L 
and a spin of S resulting in a total angular momentum of J — L-\-S 
for the electron. In the presence of a magnetic field, each state with 
an angular momentum of J splits into 2 j +1 states each with quan¬ 
tum number M. This splitting of energy states is proportional to 
Bg where g is the lande factor and depends on the energy levels, 
orbital and spin angular momentum (Reiners 2012). Transition be¬ 
tween two states must obey the selection rule AM — —1,0,1 and 
so there are three energy levels at present of magnetic field. Now 
we assume Fo( A) be a profile of a sensitive line in the absence of 
magnetic field, when a magnetic field switch on, it changes to 

F m (\) = F 0 (A)*(|«5(A —AA) + |5(A + AA) + 6«5(A)) , (10) 

where * indicates convolution and F m (A) and £(x) are spectrum 
profile in the presence of magnetic field and Dirac delta function 
respectively. In above equation a and b are two constants which 
depend on direction of magnetic field and line of sight ( 7 ), and are 
given by (Bobcock 1949): 

a = 0.5(1 + cos 2 7 ) , b = 0.5 sin 2 7 . (11) 

The wavelength displacement due to magnetic field is: 

aA = 36.8 4 (X) ( _^_)>( I ^)„A ( 12 ) 

Hare Ao is the rest frame wavelength of the line. Converting to 
velocity we have: 

A “ = 1MI S )( a7 )( li5 )Kms " <13) 

The typical magnetic field on the surface of stars is on the order 
of a kilo Gauss which results in velocity shifts of the order of km 
s - 1 for Ao ~ 0 . 8 /im. According to Eq. (13) infrared and longer 
wavelength transitions have larger shifts than those of optical. For 
example 1.5 6/im of (Fel) and 2.22 gm of (Til) are two sensitive 
line which can be used in infrared to detect a magnetic field. 

The brightness profile of a star can be written as: 

P( A) = F( A) * V(X) * M(A) * /(A) , (14) 

where F( A), V(A), M(A) and /(A) indicate the intrinsic profiles, 
velocity, micro and macro turbulence and the telescope broadening, 
respectively. In Fourier space Eq. (14) transforms to 

p(k) = f(k)v(k)m(k)i(k) , (15) 

where convolution converts to ordinary multiplication in Fourier 
space. The intrinsic profile is 

F(\) = fF m (\) + (1 — f)Fo(X) , (16) 

and using Eq. (10), it transforms to 

f(k) = / 0 (/c)[/a(cos(27rAA/c) - 1) + 1] . (17) 

Our strategy is to measure the magnetic field base on the Robinson 
method (Robinson 1980). According to this method, by dividing 
the line profiles of two transitions in Fourier space we have: 

p(gc,k) _ /(g 1 ; fc) 
p( 92 ;k) f{g 2 ;k) ’ 


MNRAS 000, 000-000 (0000) 


























4 Mehrabi et al. 


where all broadening agents from Eq. (14) will be omitted. Using 
two lines with small (insensitive) and large (sensitive) lines one 
can easily separate magnetic field signal from other broadening 
agents. Our following work on the measurement of magnetic field 
essentially base on the Eq. (18). In this equation f(g,k) comes 
from Eq. (17) and depend on the fo(k) and AA where fo(k) is the 
Fourier profile for B = 0 and A A as a function of g that is given 
by Eq. (12). To find f(g\k) we need to calculate Fourier transform 
for a discrete sample. According to Nyquist theorem if sampling 
rate in velocity space be Av, its Fourier transform is confined in 
the range[— 2 ^ 7 , 7 ^ 7 ]. This limitation imply that we can’t detect 
magnetic field with A A/A << A v/c, as we expected. 


4 DETECTABILITY OF THE MAGNETIC FIELD AT 
CAUSTIC CROSSING 

In this section we study the detectability of the magnetic field at 
caustic crossing. To do so we consider two sets of spectroscopic 
observations as following: 

(i) spectral resolution of Av = 1 Km s -1 with a signal to noise 
ration (SNR) of 20 

(ii) spectral resolution of Av = 5 Km s -1 with SNR = 50 

and perform a Monte Carlo simulation to study the detectability 
of magnetic field at caustic crossing. These are a typical value 
for a spectrogragh and in this paper we aim to study the effects 
of caustic crossing on field detection. We simulate several events 
with different parameters (q,d,uo,tE) and near central caustic 
(high magnification), the detectability of a magnetic field is cal¬ 
culated. Moreover in this study we assume that the spot can lo¬ 
cate over each part of the source with equal probability so the 
location of spot is chose randomly in each event. Now we need 
to define a criteria to select detectable events. Since the signal of 
magnetic field comes from Eq. (17), by Assuming a sensitive line 
with (Ao = 0.8 fim, g 1 = 2.5) and an insensitive line with 
(Ao — 0.8 /im, g 2 — 1 ), we have: 

/(gi;fc) _ /o ew (fc) [/a(cos(27rAA 3en fc) - 1) + 1] 
f( 92 ',k) fo ns (k) [fa(cos(2-KAXi ns k) — 1) + 1] 

where superscript sen and ins stand for sensitive and insensitive 
lines, respectively. In the following we assume that fo en (k) and 
frik) are known function that in principal can be obtained by 
solving the radiative equation of transfer (Saar 1988). In this case 
the signal comes from function 

c/ o\ = /«(cOs(27rAAsenfc) ~ 1) + 1 .... 

{ ' fa(cos(2nAX ins k) - 1) +1 ' ^ ’ 

The value of a is given by Eq. (11) and we assume an average 
value 7 = 34 from (Marcy 1981). For B = 0, we have S(B = 
0) = 1, so any deviation from 1 indicates 77 % 0. We assume a 
Gaussian distribution for flux errors which is distributed like true 
space in Fourier space, then among all simulated events, those with 
signal more than 5 a are flagged as detectable events. We should 
notice that in our analysis / is a function of time and we use a mean 
value for this quantity (see Fig. 1). This mean value is obtained by 
averaging f(t) around caustic line. 

The magnetic field detectability is defined as a number of 
events which pass our criteria per total number of events and the re¬ 
sults for strategy 1 and 2 are shown in Figs. (3) and (4). In strategy 
1 , for / = 0.001 the detectable events are very rare and for very 
high magnification the detectability is less than 20%. For a spot 
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Figure 3. Percent of detectable event in the A max — B plane. Dark area 
(yellow point) present undetectable (fully detectable) events. The size of 
spots is presented above each diagram. 
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with initial size / = 0.005 magnification larger than around 500 
is required to detect magnetic field but for / = 0.05 events with 
small magnification are fully detectable due to high signal from a 
large spot. The results for strategy 2 are almost like the results of 
strategy 1 but due to low spectral resolution, small magnetic field 
can’t be detected while high SNR help to detect magnetic field for 
/ = 0.001 in contrast with strategy 1. 

Since in a microlensing event it is possible to realize type of a 
source star before caustic crossing, selecting a more sensitive star 
leads to a higher chance to detect magnetic field. M-type stars have 
radii smaller than solar radius and also rotate slower so these condi¬ 
tions lead to produce a large field in such stars. In addition two un¬ 
wanted broadening mechanism, rotational and temperature broad¬ 
ening, in a M-type star are smaller compare to a sun like star. All 
of these effects will facilitate field detection in an M-type source 
star, (for more information on the magnetic field of M-type star see 
(Reiners & Basri 2007, 2010)). Sun like stars may be active or less 
active but it is not possible to find out that a sun like source is active 
before caustic crossing, (to see a list of field detection on a sun like 
star refer to (Reiners 2012) ). Moreover giants stars can be rather 
active. The field of giants are much weaker than dwarfs but the size 
of spots may be very large (Petit et al. 2004). Giants stars also are 
good candidate to detect magnetic field. 


5 CONCLUSION 

Several teams including OGLE, EROS, MOA, PLANET and Mi- 
croFUN are observing sky to find a microlensing events. These sur¬ 
veys detect several thousand events per year and a fraction of them 
are binary system. To our knowledge there was no attempt to study 
the magnetic field of source star in these events. In this work we 
propose an interesting approach to detect magnetic field which may 
be existed on the surface of the source. Our detection mechanism 
based on the Robinson method which investigate two line profiles 
of a star in Fourier space to deconvolve all broadening agents from 
magnetic one. We study the possibility of magnetic field detection 
when a source is crossing the caustic line. At a caustic crossing 
not only the flux of a source increase which results higher SNR, 
but also the coverage fraction of spot increases due to magnifica¬ 
tion gradient. We consider two spectroscopic strategies which are 
in agreement with recent powerful spectrograph, and find the de¬ 
tectability of magnetic field in the range B £ (0 — 5)kG and 
magnification in the range /i max G (10 — 2000). In addition, in 
a microlensing events, is is possible to find out the type of source 
stars before caustic crossing. This situation provide possibility of a 
more active source selection. Among all stars, M-type stars and gi¬ 
ants are two good candidate to detect a magnetic field. Our results 
show that high magnification channel provide an unique opportu¬ 
nity to study the magnetic field of stars which are located at bulge 
and disk. 
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